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ANNOTAi'ION 


An analysis of kinetic equations which describe, in par- 
ticular, the establishment of the spectra of Langmuir turbu- 
lence is presented. For the analytical nuclei and linear in- 
crement, a stationary spectrum. If it exists, is concentrated 
on a variety of dimensionality, the smallest dimensionality, of the 
initial problem at. points, on closed lines or surfaces. In 
the case where stationary distribution consists of many peaks, 
with a sharp inclusion of the increment of pumping in the sys- 
tem, secondary turbulence occurs. The position of peaks is 
established but their amplitudes complete undamped (in the ab- 
sence of noise) oscillations. It is pointed out that establish- 
ing spectra can occur only during adiabatic inclusion of pump- 
ing. Then, one realizes a quasi stationary condition. It is 
significant that the adiabatic condition here is more rigid 
than the ordinary (Tt>>l, y — the rate of pumping of oscil- 
lations, t — the time of inclusion of instability) by several 
hundred times. 
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ESTABLISHMENT OP THE SPECTRA OP KINETIC TURBULENCE 

by 

L. A. Bol'shov, A. M. Dykhne, V. P. Kiselev 
and A. Kh. Pergament 


1 . Introduction /5 * 

In a number of physical objects described In an approxi- 
mation of a solid medium, with a fairly high level of excita- 
tion, complex motion occurs which is described in the language 
of mutually interacting degrees of wave freedom; the Langmuir turbu- 
lence in plasma, acoustic turbulence, noise in amplified media, spir 
waves in magnets, etc. Motion of this type, with not too great 
intensity of turbulent noise, can be described approximately ! 

as weak turbulence, as a kinetic equation for interacting waves j 

with random phases. Interaction of the waves with nondecay by i 

f 

the principle of dispersion is accomplished in the lowest mag- | 

nltude according to Intensity of the wave by induced scatter- | 

ing on other exlcitations (in the case of Langmuir turbulence 5 

in an Isothermal plasma, scattering of the waves occurs basic- | 

ally In ions). i 

A kinetic equation for the number of occupations n(k,t) I 

is well known [1,2] 

«!. »2 ( 1 ) 

here is an incremem, (or decrement) of surging of waves in t 

a linear approximation, Aj^j^i is a nucleus, which describes the I 

process of the scattering; is the intensity of the source 


^Numbers in the margin Indicate the pagination in the foreign 
text . 


caused by thermal noise and Induction processes of a higher i 

order. 

1 

i 

In situations which are of interest from a point of view j 

t 

of heating plasma with an electric beam, or by electromagnetic \ 

radiation, the intensity of excited waves is much larger than } 

the level of thermal noise. The latter added to (1) then is | 

considerably smaller than the other members. There is inter- '! 

est in the study of the spectra of turbulence in time periods smaller ! 
than The behavior of spectra for large time periods is con- 

.Tidered in references [3,6]. ^ 

The question of establishing the spectrum of Langmuir 
turbulence of plasma was studied within the framework of 
equation (1) in a number of works [3-6]. How'ever, the study 
was Incomplete [3], or contained a number of erroneous con- 
clusions [5]. 

We wish to approach the study of the spectrum from both 
positions, with the possibility of not using the specifics iof 
Langmuir turbulence, expressed in the actual form of the func- 
tlorBYi(, and . We will try to answer the following ques- 
tions : 

1) What are the possible types of stationary spectra? 

2 } What are the characteristic traits of the spectra 

in the last stages of evolution? 

3) Which can be nonstationary spectra and what are the 
causes for their establishment, or nonestablishment? 


2. Stationary Spectra . 

Stationary spectra, if such exist, satisfy the equation 




(2) 


Here the nonlinear Increment 

'V • arn ♦ 

It is possible to solve equation (2) of three types. 

(a) with all ^ (solid spectrum). 

(b) in certain fields 5 — space. 

(c) for multiple small dimensions (on surfaces, 
on curves and at separate points ^ — space). We will con- 
sider these possibilities in more detail. 

(a) Although examples exist in which there is the pos- /7 
siblllty (a) (for example [5]), this possibility must be 
designated as exclusive. Moreover, in this case when all ^ 

it must satisfy the Integral equation The so- 

lution of this equation, generally speaking, has a variable 
sign. The natural physical requirement for the positive 
state of n^ can be satisfied only with a special coordinated 
selection of the function and . 

(b) In this case, there must be rj^=0 in fields where 
the spectrum is concentrated (that is, Hj^/O) and Tj^<0 (for 
stability of a given spectrum) in all other regions of ic- 
space. We note that the function which satisfies this con- 
dition Is an insignificant function of its independent vari- 
able It is easy to see that any non-analytlclty and even 
break in the function n(Jc') introduced into Fj^ cannot result 

in the non-analytlclty of Fj^ (inasmuch as n(i<» ) is Integrated). 
Thus, the single cause of non-analytlclty of Fj^ can be non- 
analyticity of the nucleus or the linear Increment Yj^. 

The ’’striped" spectrum can be stabilized if or yj^ have 

characteristics as functions of their independent variables. 



A model example with a nucleus can Illustrate such a 
spectrum 


* 


f i . loa' 

' ^’when 


-i when 


( 4 ) 


If y(k) has one maximum at point and Y(kQ)®-Y(+*) > then 

a spectrum is realized oncentrated on semistraight 

fl s J 'f ( 5 ) ^ 

* [ 0 vflien M<M0 


The peculiarity in nj^ when k«k^ involves the peculiarity of 
the nucleus. 


It is interesting to note that for a nucleus which de- 
scribes the scattering of the Langmuir waves In ions, with 
the use of which the Langmuir turbulence was studied in 
references [3-6,8], the nucleus can be both analytical 

(if the spectrum is concentrated in the k>>kp/m7M field, in 
this case scattering is. significant at large angles), and 
nonanalytlcal (if the spectrum is localized in the field of 
k<<kj^/m/M, where small angle scattering of the plasmons plays 
a significant role). 


(c) In the case of analytical functions and Yj^ the 

most typical is the third variation in which the spectrum is 

concentrated to a multiplicity of dimensionality, the smallest dlmen 
sionality of the initial problem at points, on lines or sur- 
faces. The existing lines or surfaces inevitably must be 
closed because in the opposite case the peculiarity of is 
at the terminal points (lines) of the corresponding multiplicity . 

H 


Thus, the spectrum cannot be concentrated on sections of the 
beams (notwithstanding the statement contained in reference 
[ 5 ]). In the variation of the function considered, on the 
aforementioned manifolds it has a maximum equal to zero and 
negative in all other space. In reference [53 the regular or 
singular character of the spectrum Involves the width of the 
increment In relation to the width of the nucleus. As 

is apparent from what has been presented above, the character 
of the spectrum is defined exclusively by the analytical 
properties of the Aj^j^, nucleus and the linear Increment of y^. 

During an analysis of establishing the spectrum of type ^ 
(c) in a case of Langmuir turbulence, it is necessary to keep 
in mind that with adequate constriction of the peaks they can 
become nonstationary as to modulation [ 73 . This circumstance 
lindts the use of the picture considered here to large time 
periods . 

Later on we will consider only the (c) type spectrum. 

When studying later stages of evolution of the spectra of this 
type, for simplicity, we will limit ourselves to a unidimen- 
sional case where the spectrum will be localized at N points 
x^(i=l,. . .N) 

♦ •f 

Then only when x=x^ (i=l , . . . N) r(x)=0 and from requirements 
for stability of the given spectrum r'(Xj^)=0, and r(x)<0 at 
all remaining points. For determining a given spectrum it 
is adequate to solve the system: 

I jV*i)+S » 0 

»>i 
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The peculiarity of system (7) involves the fact that not 
only the position of the peaks of and their amplitude 
are unknown but also the number of peaks N [illegible]. With 
weak instability (when the 7 increment is small and positive 
only in a narrow field according to x) the distribution estab- 
lished consists only of one peak. As the Increment Increases 
the stability condition breaks down (r(x )<0 everywhere except 
at one point) and at this point where this occurs, a second 
peak of the established distribution occurs. With a further 
increase in y» one can find a third peak or even more. 

A numerical determination was made of this stationary /lO 
distribution for the nucleus 

${oc) • X tacp/-oc*/2) ( 8 ) 

(this nucleus makes the transition in the direction of large 
x) and two model Increments y(x) 

1) ir^x) 

2 ) i>i 

stationary distribution consists of one peak when 
0<a<2.21 and 1<B<1.314. When 2.2l£a<3-3 and 1 . 31^lB<l . 72^ , 
the established spectrum contains two peaks. When a=3.3 and 
B-1.72A, a third peak is apparent and so forth. The coordin- 
ates and amplitudes of the first three peaks in relationship 
to a and B are presented in Figures 1-2. 

3. Evolution of the Structure of a Separate Peak . 

Now we will study the process of the tendency toward 
stationary distribution of type (c). Here one can Isolate 
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/23 



a 



b 


Figure 2. Relationships of the amplitude /20 

of peaks of stationary dlstri- 
bution to the level of pumping 
with model increments. 

a) Y (x)«x(a-x) 

b) Y (x)-1/(1+x2)-1/B 


two questions. The firbt of these questions is on the struc- 
ture of a separate peak in the later stages of evolution and 
the second question la on the behavior of the spectrum as a 
whole. When studying the structure of a separate peak, for 
simplicity, we will limit presentation to a case where sta- 
tionary distribution contains a total of one peak. This situ- 
ation occurs with a fairly weak Instability when the increment 
Y (x) is small and positive in the narrow field of x. In this 
case, the stationary condition consists of one peak located at 
a point where y(x)»0. For large time values, the distribution 
of n(x) which tends toward a 5-form, becomes fairly narrow. /II 
This makes it possible to expand the G(x-x’) nucleus to an Integral 
and also the increment y(x) Is close to the stationary position 
of the peak. For the Gaussian nucleus 


• X - oc*/z 

y(x) « 96 


( 10 ) 


( 11 ) 


The solution of equation (1) is written in the following form 




( 12 ) 


with calculation of (11), (12), equation (1) can be rewritten 

Multiplying (13) by 1, (x-xi), (x-Xi)2 and Integrating by x 
we obtain 


V, • ^ .4 (x!-* 
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( 15 ) 


^ - t/lac, A* H ♦ 

A • -^A* 


The established distribution corresponds to 


V«*< 


*ss 


Prom (16) we have 


A*- 


( 16 ) 


(17) 

( 18 ) 


Considering that and x differ little from stationary 
values, let us assume Nj-N^-N^g and leaving only linear mem- 
bers in (1^) and (15) along Nj and Xj we have” with a calcula- 
tion of (18) 

( 19 ) 

i • AC« 

( 20 ) 

By introducing z-xjt (21) from (19), (20), we find /12 


from which 


that is 


“ v/»e “ ^ 

(22) 


(23) 


(2iJ) 

-/'•S ^co^( ) 

(25) 


Thus, the amplitude of the peak and the position of its 
center contain oscillating additions which decay with time. 
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The oscillations of the position of the center of the packet 
decay most rapidly. The width of the packet decays 
like t“^/2 and finally oscillation of the height of the peak 
(<\.t**^/‘*) decays most slowly. 

A numerical study was made of the process of making a 
permanent establishment in this case. Equation (1) has the 
form 

where y(x)=x(1-x). The level of thermal noise was selected 
as the initial level of n(x). (It appeared that the behavior 
of solving equation (26) hardly depends on the initial level 
of n(x).) Calculations were made with different noise levels 
(during observation of the condition £j^<<max(Y) ) • We at- 
tempted to show the characteristics of the establishment 
process not depending on the noise level. 


With numerical integration according to time of the kinetic 
equation ( 26 ) the known explicit system with a second magnitude 
of precision was used which for equation 




(27) 


can be writeen in the following form /13 


The integral in the right section of (26) is calculated accord- 
ing to the Simpson method and then preliminarily calculated 
values of the nucleus for different values of x are used. We 
note that n(x) was calculated as changing in time only in a 
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certain limited field of values of x within which the full 
increment r(x) is deliberately negative for any values of 
t and where one assumes n(x)=ej^. 

The process of showing the solution for a single peak 
permanent establishment has the following character. At 
first, in the field of the maximum increment, a broad peak 
develops (on the order of the width of the increment) which, 
having reached a fairly large value (approximately lOOvmax) 
begins to shift toward the large values of x (by virtue of 
selecting the sign of the nucleus). Having overshot the 
position of permanent establishment, the peak is damped and 
shifts toward small values of x, and begins to increase once 
more when x is smaller than the stationary value but larger than 
X corresponding to the maximum increment. Evolution n(x) in 
time is presented in Figures 3-^. Establishment of the posi- 
tion of the peak occurs most rapidly (when Ej^= 0.01 for yt'^lOO, 
when Ej^=0.001 for Yt'''200). The character of the appearance 
of the position of the center of the packet on a permanent 
establishment is shown in Figure 5 . After establishing the 
position of the center of the packet, decay of the peak occurs 
with a simultaneous increase in its height (see Figure 7). 
Establishing the area of the packet occurs more slowly. With 
precision less than 1?, the area S=/n(x)dx appears on the-* 
permanent establishment when Ej^= 0.05 for time Yt'''200, and 
when Ej^=0.01 for yt'''600. Thus, the rate of appearance of the 
area of the packet on the permanent establishment strongly de- 
pends on the noise level. The half-width of a distribution / 
n(x), calculated according to the formula 

(29) 

ivf nnv 



Figure 3- Evolution of n(x) for small /20 

time values . The curves 
correspond to times: 1-Yt=15, 

2-Yt=30, 3-yt=45, 4-Yt=60, 

5-Yt=75. 

hardly depends on the level of noise and beginning with /IH 

Yt'b50 is subject to the law f-i/2. (Figure 6 shows the re- 
lationship of in 1/a to s.nt.) 


One should note that the distribution of n(x) obtained 
numerically is essentially non-Gaussian. It is fairly sym- 
metrical relative to the center, but its excess SI 

Sis — (30) 


appears to be several times larger than the excess for Gaussian 
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Figure 4 . Evolution of n(x) for large /20 


time values . The curves 
correspond to times: 1-Yt=100, 

2-Yt=200, 3-Yt=300, ^1-Yt=400, 

5-Yt=500 . 

distribution equal to three. Thus, the asymptotics of (I8), /l 4 

( 2 ^), (25) catch only a qualitative character of the solution 
on permanent establishment. 

It is significant that the area of the peak like the 
half-width is established slowly which, strictly speaking, 
causes doubt as to the existence in this case of the first 
stage of the process in which the energy balance is estab- 
lished between pumping and decay according to reference [ 8 ]. 
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Figure 7. The typical phase curve of 
establishing a two-peak 
permanent establishment 
with abrupt inclusion of 
the increment. 


Setting up a Multipeak Permanent Establishment 
We will study the character of establishment in this case 
where the distribution established consists of several peaks. 

If one assumes the existence of permanently established distri- 
bution, then one can describe nonstationary evolution of distri- 
bution for fairly large time periods. For this, we will assume 
that distribution of density can be prevented as a superposi- 
tion equally depending on the time section <nj^(t)> (which tends 
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toward stationary distribution) and the oscillation addition /15 
f»j^(t). Prom formulas (2), (3) we find 

<«.«» (31) 

The first three multipliers in the right part of expres- 
sion (31) for large time intervals are a stationary distribu- 
tion which, in accordance with what has been said above, con- 
sists of a series of 6-shaped peaks. Inasmuch as the fourth 
multiplier in (31) is limited in time (due to the oscillating 
character of fij^) we have when t-*-» 

'6 

HttC'i) -* (32) 

If a stationary solution exists, then even with the presence 
of an oscillation in the amplitudes of the peaks averaged for 
time, the value of nj^ coincides with nst^ which automatically 
is completed if 

( I • i (33) 

For studying the general picture of establishing the 
spectrum for large time intervals, one can use the peak 
character n^^(t). Because the position of the peaks is 
established most rapidly, we will assume that nj^(t) consists 
of N oscillating peaks with amplitudes Na, found at points 
corresponding to a stationary distribution. In this case, 
equation (1) is transformed to a system of nonlinear equa- 
tions, the first of which Is presented in reference [5] 

( 31 ) 

Here is the amplitude of a peak with the number a. In 
the case of a random number of peaks, the character of motion 
of the system close to the permanent establishment depends on 
the number of nondamped oscillations. To answer the question 
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as to the number of nondamped oscillations in a general case, /16 
one mainly uses a linear theory. However, with a fairly 
large number of peaks, this Is difficult to do in actual 
practice. An exception is the case of interaction of waves 
which describe an antisymmetrical matrix A. From the fact 
that SpA«0 and ReA.<0 (a condition of stability of a station- 
ary distribution, — the intrinsic number of A) substan- 
tiation follows for all Intrinsic frequencies of the system 
(absence of damping in all modes). These intrinsic nondaunp- 
ing modes can be excited in the system to a different level. 

For large numbers of such peaks this phenomenon is a singular 
secondary turbulence. It can be characterized by the func- 
tions of distribution according to the degree of excitation 
of these new modes. With a fairly high degree of excitation, 
the modes Interact strongly and must be used by a system of 
nonlinear equations (31). The unusual character of this type 
of turbulence involves the fact that in a strongly perturbed 
distribution system, motion is realized with a finite number 
of degrees of freedom (on the order of the number of peaks). 


The question as to excitation of secondary turbulence is 
very complex and requires a separate study. Here we will limit 
ourselves to separate qualitative expressions. Initial in- 
stability leading to primary turbulence (positive nature of 
Y at certain k) at the same time is included with a certain 
finite velocity. We will consider two separate cases: in 

the first of these the time of inclusion of the increment 
significantly Increases the time for development of primary 
turbulence. 


jL -i- 


(35) 


In this case at any value of y» the system succeeds in 
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becoming stationary and secondary turbulence does not develop. 
In the opposite extreme case of sudden Inclusion of Instab- 
ility 

(36) 

strongly developing secondary turbulence must occur. 


We will study the simplest case where the matrix A 20 In 
(3^) Is antlsymmetrlcal and has a total of two peaks. Then 
we have 


(37) 

For this system 

Hl , JliSikzad 

Jit, Nt(Ht-Hu) (38) 


that Is, 


iit Mg ~ Mgg ittMgm 


(39) 


The phase trajectories of the (37) system are closed lines with 
the center at a stationary posltlonNig, N 25 . Thus, the system 
completes non-damping oscillation, not tending toward the 
permanent establishment. The amplitude of these oscillations 
depends on the Initial conditions; then In the case of zero 
Initial conditions, the span of oscillations Is great. It Is 
easy to prove that In this case the mean value of the exponent 
In (32) actually equals one. This condition is neutrally 
stable — small excitation leads to a system with one cycle 
close to another. The Introduction of noise Into the right 
part of the model system of ( 37 ) results In the appearance of 
a focus Instead of a cycle and to the establishment of 
equilibrium. The rate of approach to a permanent establishment 


Is Increased with an increase In noise. 


A numerical study was made of the process of establishing 
a solution for equation (26) when 

-i// (HO) 

which when b*»1.5 corresponds to stationary distribution from /l8 
two peaks. Also the time of Inclusion of the increment of 
instability was varied, producing 

With an abrupt Inclusion of the increment (tg»30) at first 
fairly rapidly (when Yt'v60) at points corresponding to sta- 
tionary distribution, two fairly broad peaks are established 
whose area varies and whose width gradually contracts. It 
is convenient to draw the phase trajectories of such a process, 
having applied along the axis of the coordinates the area of 
the first and second peaks. A typical example of this phase 
trajectory is shown in Figure 7. The system is completed by 
slowly damping oscillations around a stationary position. 

With slow inclusion of the increment (tg»200 and 
especially tg^^OO) the character of the process of establish- 
ment sharply changes and a quasi stationary condition for 
establishment is realized. Figure 8 shows phase trajectories 
of a system for different time inclusions (tg^SO, 200, ^00, 
respectively). Numerical calculations show that in inequalities 
(35), ( 36 ) it is necessary to substitute one for the small 
parameter 'v-lO”^. This apparently involves the fact that for 
fairly narrow peaks n(x) differs noticeably from zero only at 
points where there is a small full increment r(x) (because 
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Vf;, 


|>V'^ 






\, V 




Pig. 8a 


Pig. 8b 


Fig. 8c 


Figure 8. The relationship of the character 

of establishing a two-peak permanent 
establishment of the time of Inclu- 
sion of the Increment, time of In- 
clusion: a) to«30, b) tn»200, 

c) tn*^00 







In the permanent establishment r(x)«o). 

It is natural that the rate of convergence toward a 
permanent establishment with an abrupt inclusion of the in- 
crement is increased with an Increase in the noise level. 

Figure 9 shows phase trajectories of the system when Cj^*0.01 
and Cj^"0.001 (tg"30). Then, for a low noise level (cj{*0.001) 
the sections of the phase trajectory when yt>200 are close to 
the maximum cycles of the model system (37) with the correspond- 
ing initial conditions. 

5. Conclusions /19 

1. The character of the established spectrum of kinetic 
equation (1) is determined by the analytical properties of the 
nucleus and the linear increment Yj^. In a case of analytical 

\k’ ^k* stationary spectrum. If it exists, is concen- 

trated on varied dimensionality, the smallest dimensionality of 
the initial problem, at points on closed lines or surfaces. 

2 As an example of establishing a single peak permanent 
establ*’/iment analytically and numerically it is pointed out 
that the height of the peak and the position of its center con- 
tain oscillation additives which decay with time. At first, 
the position of the peak is established significantly smaller 
than its area. Then distribution of the intensity with'n the 
peak is essentially not Gaussian. The time of establishing a 
permanent establishment is several f undr .o of times -’ger than 

^/“Ymax* 


3. In a case when stationary distribution consists of 
many peaks, with the abrupt inclusion of the pumping Increment 
into the system, secondary turbulence occurs. The position of 
the peaks is established and their amplitudes complete non- 
damped oscillations. In the distribution system considered. 










a finite number of degrees of freedom is excited (on the 
order of the number of peaks). 

As an example of establishing two-peak stationary 
distribution, the question of the character of secondary 
turbulence Is studied in detail. It is pointed out that with 
abrupt Inclusion of the increment, nonstability of the system 
completes non-damped (in the absence of noise) oscillations 
with large amplitude around a stationary position. During 
adiabatic inclusion of the Increment (when the condition of 
adiabatlclty is very rigid ^ ^<10”*) a quasi stationary 
state of establishment Is realized. 


2^1 


